Interaction effects and quantum phase transitions in topological insulators 
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We study strong correlation effects in topological insulators via the Lanczos algorithm, which we 
utilize to calculate the exact many-particle ground-state wave function and its topological proper- 
ties. We analyze the simple, non-interacting Haldane model on a honeycomb lattice with known 
topological properties and demonstrate that these properties are already evident in small clusters. 
Next, we consider interacting fermions by introducing repulsive nearest-neighbor interactions. A 
first-order quantum phase transition was discovered at finite interaction strength between the topo- 
logical band insulator and a topologically trivial Mott insulating phase by use of the fidelity metric 
and the charge-density- wave structure factor. We construct the phase diagram at T = as a 
function of the interaction strength and the complex phase for the next-nearest-neighbor hoppings. 
Finally, we consider the Haldane model with interacting hard-core bosons, where no evidence for a 
topological phase is observed. An important general conclusion of our work is that despite the in- 
trinsic non-locality of topological phases their key topological properties manifest themselves already 
in small systems and therefore can be studied numerically via exact diagonalization and observed 
experimentally, e.g., with trapped ions and cold atoms in optical lattices. 
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I. INTRODUCTION 

The conventional definition of a topological insula- 
tor (TI) relies on an analysis of non-interacting band 
structures and from the mathematical point of view 
represents topological classification of the spectrum of 
matrix Hamiltonians. 1-7 While this approach is non- 
generalizable to interacting systems, it is clear that they 
may exhibit physical phenomena associated with non- 
trivial topological spectra. Qi and co-workers recently 
suggested a more general definition of a topological insu- 
lator, potentially suitable to interacting systems, which is 
defined as a state where charged carriers give rise to axion 
electrodynamics with a non-trivial electro-magnetic cou- 
pling term. 7,8 One possible approach would be to use the 
Volovik formula for topological indices, which involves 
exact Green's functions. 9 In addition, a technique uti- 
lizing a Green's functions formula for time-reversal in- 
variant topological insulators was recently presented. 10 
Finally, the topological properties can also be studied 
via explicit analysis of the edge states and the density 
of states (DOS) spectrum with an eye on the localized 
boundary modes that would connect bands across a gap. 

Interacting topological models are expected to be 
even richer than the non-interacting systems and may 
host a variety of different phases and phase transi- 
tions. However, despite this expected rich variety of 
phenomena and fundamental interest, the effect of in- 
teractions in topological insulators has remained largely 
unexplored and the current knowledge of this issue is 
rather limited. u ~ 16 The existing open questions include 
the following, (i) Can topological properties be tuned 
by interactions? (ii) What is the nature of the quan- 



tum phase transitions separating phases from different 
topological classes? (iii) Arc topological Mott insula- 
tors possible, where the non-trivial phase would arise 
entirely due to interactions? And an even more provoca- 
tive question: (iv) can topological Mott insulators exist 
in bosonic systems? The main difficulty in addressing 
these questions is that they all involve strong correlation 
physics, where essentially no quantitatively reliable ana- 
lytical methods exist. The usefulness of quantum Monte 
Carlo techniques 17 is expected to be limited by the infa- 
mous sign problem. 18,19 The last resort resides in exact 
diagonalization methods and other unbiased numerical 
techniques, which arc limited to small system sizes as 
the dimension of the Hilbert space grows exponentially 
with system size. The latter circumstance is important, 
particularly because topological phases involve non-local 
effects and it is not clear if they survive in small systems 
[e.g., due to detrimental effect of tunneling between the 
edge states in finite lattices with open boundary condi- 
tions (OBC)]. 

In this paper, we employ the Lanczos algorithm 20,21 
to calculate the exact ground-state wave function of in- 
teracting lattice models which are known or expected to 
have a topological ground state. We focus our analysis 
on four observables that provide unique insights into the 
properties of these lattice models. The topological phase 
is examined with the local density of states (LDOS), 
which we determined for both open and periodic bound- 
ary conditions (PBC). As we will discuss in detail be- 
low, the key topological features are still discernible in 
relatively small clusters (at least in two-dimensional lat- 
tices), which are within the reach of exact diagonalization 
techniques. The next two observables of interest are the 




FIG. 1: (Color online) Schematic V-<f> phase diagrams for 
spinless fermions in the Haldane model at half filling. For 
< (f> < T) the system is a topological insulator (TI) at weak 
coupling and becomes a Mott insulator at finite interaction 
strengths via a first-order phase transition. At cj> = and <j> — 
n, the TI phase is not possible. Instead, the phase at weak 
coupling is a semimetal. The points in the phase boundary 
were calculated for a 24-site cluster (24D). 



edge currents and the structure factor, which we define 
carefully in the next section. The final observable of in- 
terest is the ground-state fidelity metric g, defined below 
in Eq. (10). The fidelity metric is related to the rate 
of change in the overlap between the ground-state wave 
function in Hamiltonians that differ by a small control pa- 
rameter. The ground-state fidelity was originally studied 
in quantum information theory and has been shown to be 
a sensitive indicator of quantum phase transitions. 
In conjunction with the structure factor, the fidelity met- 
ric is used to characterize the nature of the transition 
from the topological insulator to the topologically trivial 
Mott insulator. 

To test the feasibility of the Lanczos algorithm and to 
explore the main questions, we chose the Haldane model 
on a honeycomb lattice, but with-nearest-neighbor re- 
pulsive interactions. The non-interacting base model is 
known to have a simple topological phase transition to 
a trivial band insulator, 1 which can be tuned by adjust- 
ing a staggered chemical potential, M. Moreover, in the 
limit of weak interaction strengths, the bulk gap and 
non-trivial topology are known to persist. 32 ~ 34 This al- 
lows us to calibrate the techniques on the non-interacting 
and weakly interacting fermion model, whose topological 
ground state is known and well understood. The main 
conclusion of this initial analysis is that the local density 
of states for the bulk and the edge are the most useful 
metrics to probe the topological nature of the ground 
state. On the other hand, we find that the existence of 
the chiral edge currents may be misleading, as those may 
be present even in topologically trivial phases and their 
origin and stability cannot be determined just by looking 
at the current texture in the ground state. Indeed, this 
result has been independently verified for non-interacting 
systems. 35 

With these caveats in mind, we analyze the strongly 
interacting Haldane model with spinless fermions. The 
model is characterized by four independent parameters: 
nearest-neighbor hopping amplitude, t l5 next- nearest- 
neighbor hopping amplitude, £2, and its complex phase 
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FIG. 2: (Color online) Illustration of (a) the Haldane model 
and (b) the modified Haldane Hamiltonian. In (b), the direc- 
tion of positive complex phase in the next-nearest-neighbor 
hopping on the sublattice denoted by the squares is reversed. 



4>, and the nearest-neighbor coupling, V. We limit our 
analysis of the interacting model to the case of zero stag- 
gered potential, M = 0, which apart from the trivial 
case of <j) — (which is a semimetal) is guaranteed 
to be a topological insulator for V <C £1,2- We find 
that as V is increased, the system undergoes a topolog- 
ical quantum phase transition into a topologically triv- 
ial Mott insulator with the critical interaction strength 
V c ((f>) = V c (ir — <f>) monotonically increasing increasing 
from V c (tt/180) = 0.19ti to V c (tt/2) = 2.29£i within half 
a period (see Fig. 1). The appearance of the Mott insu- 
lator is not surprising, but it is not a priori clear if the 
Mott transition must always coincide with the loss of the 
topological order or if there are two separate transitions. 
Below, we show compelling evidence that the former sce- 
nario is the one realized in this model, that there is a 
single first-order quantum phase transition at intermedi- 
ate couplings. 

Finally, we populated the Haldane honeycomb lattice 
with interacting hard-core bosons in the hope that a 
phase with non-trivial topological properties may emerge 
in this Bose-Hubbard model at intermediate couplings. 
However, it was found not to be the case. At half fill- 
ing, the system remains a trivial supcrfiuid until con- 
siderable values of nearest-neighbor repulsion parameter 
is reached, at which point a superffuid-to-Mott-insulator 
transition takes place. However, it is not a topological 
Mott insulator because the DOS shows no evidence of 
an edge state crossing the Hubbard gap. Because this 
hard-core boson model has interesting features beyond 
the focus of this paper, we defer a full discussion of its 
properties for future work. 36 

The remainder of this paper is structured as follows. 
In Sec. II we define the model and the key measurements 
used to characterize the ground state: local density of 
states, structure factor, and the fidelity metric. In ad- 
dition, the Lanczos algorithm is briefly discussed. In 
Sec. Ill, we examine the exactly solvable non-interacting 
Haldane model and show that topological phases can be 
observed on small lattices. Next, we present the results 
for interacting systems in Sec. IV, with Sec. IV A devoted 
to spinless fermions and Sec. IV B to hard-core bosons. 
The main results and outlook arc summarized in Sec. V. 



II. MODEL AND MEASUREMENTS 

The Haldane Hamiltonian 1 is a well-known model of 
free fermions that features the anomalous quantum Hall 
state. In real-space, the Hamiltonian is given by 
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where c\ (c^ ) represent the fermion creation (annihila- 
tion) operators at site i and ?i; = c[c i is the corre- 
sponding number operator. Here M is a staggered po- 
tential that breaks the symmetry between the two sub- 
lattices of a honeycomb lattice, indicated by odd or 
even values of er(i), and t\ (£2) are the nearest-neighbor 
(next-nearest-ncighbor) hopping amplitudes. The next- 
nearest-neighbor hopping term has a complex phase 
4>ij = ±cj). In the Haldane model, the origin of this term 
is a varying magnetic field that has zero net flux in a given 
hexagon. This is illustrated in Fig. 2(a), where the sign 
of the complex phase is positive in the direction of the 
vector between next-nearest-ncighbor sites. In Fig. 2(b), 
we show a slightly modified version of the Hamiltonian, 
which differs from the Haldane Hamiltonian in that the 
direction of positive complex phase on one sublattice is 
reversed. 

In this paper, we investigate the properties of the Hal- 
dane model for interacting particles (spinless fermions 
or hard-core bosons) at half filling. The Hamiltonian 
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and V describes a repulsive nearest-neighbor interaction. 
To study the interacting system, we utilize the Lanc- 
zos algorithm. 20 (Details of the algorithm in the con- 
text of Hcrmitian matrices can be found in Ref. 21.) 
The main advantage of this technique is that it provides 
extremely accurate information about the ground-state 



wave function for interacting quantum systems. Unfor- 
tunately, memory limitations restrict the size of the clus- 
ters that can be studied. For the Hamiltonians of in- 
terest in this paper, we consider various clusters with 
periodic boundary conditions. The largest cluster stud- 
ied has 30 sites, which at half filling has a Hilbert space 
of 30!/(15!) 2 ~ 1.55 x 10 8 , near the limit of what is 
accessible on present-day computers. We also consider 
one cluster with 24 sites and open boundary conditions, 
which is used to study surface related effects. As dis- 
cussed above in Sec. I, the ground-state wave function 
determines the topological nature of a quantum system. 
More importantly, as we will show in Sec. Ill, the topo- 
logical properties in the anomalous quantum Hall system 
can be observed for small systems with only 24 sites (or 
even smaller) . These important features of the topologi- 
cal insulator make the Lanczos algorithm a suitable tool 
in the study of the topological properties. 

As a vital feature of the Haldane model, there is no 
net flux in each unit cell of the lattice. The absence of 
net flux turns out to be crucial in the study of the topo- 
logical edge states. As we implement the open bound- 
ary condition, this could introduce a local flux near the 
boundary. This local flux induces Faraday currents at the 
edge, which hybridize with the topological edge currents 
and prevent a clear observation of the topological edge 
states. Thus, it is crucial for the edge configuration to 
be chosen carefully to avoid the net flux at the edge. In 
the numerical works reported here, we removed the next- 
nearest-neighbor hopping terms between the sites on the 
edge that do not pass through a hexagon to preserve this 
property. 

One of the distinguishing characteristics of a topo- 
logical insulator is the presence of a (topologically pro- 
tected) conducting edge state with an insulating bulk. 
The LDOS can be measured via scanning tunneling mi- 
croscopy in conventional condensed matter systems. 
For a nonintcracting system, it is defined as 
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where | (i\^ n ) | 2 is the weight of the state |\I/„) with a 
particle in site i. In the Lanczos algorithm, the local 
density of states is given by 
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where |^„) = Yl m a m I^W * s ^ ne n ^ n eigenvector with particle number N and energy eigenvalue E^ and 
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are the orthonormalized vectors determined in the Lanc- 
zos procedure. In addition, it can easily be shown 38 that 
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where N' = N ± 1. Further details on dynamical prop- 
erties in Lanczos can be found in Refs. 38-40. 

To better understand the edge states in our models, 
we considered the current on every bond in the lattice. 
Between sites m and n, the magnitude of the density 
current J mn is defined to be 41 
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where q is the charge carried by a particle and r mn is the 
magnitude of the position vector between the two sites. 
For large repulsive interaction strengths, the topologi- 
cal insulator will give way to a trivial charge-density- wave 
(CDW) insulator through a topological phase transition. 
In the limit V — > oo, the ground state will be a perfect 
CDW, where one of the two sublattices is occupied while 
the other is empty, leaving lattice translational symmetry 
intact but breaking reflection symmetry 42 The correla- 
tion function that describes the CDW phase is 

C{v i -v j ) = {{nt-n b i ){n a j -n])), (7) 

where n° and n\ are the number operators on sublattice a 
and b in the ith unit cell, respectively. The corresponding 
structure factor is 
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Because the long-range order is diagonal, S(k) will be 
maximal at k = and we define S'cdw = S(k = 0). 

In order to study this quantum phase transition, we 
introduce an observable related to the fidelity F. Let 
|*o(V)) be the ground state of H(V) and \V (V + SV)) 
be the ground state of H(V + SV). The fidelity F(V, SV) 
is simply the overlap between these two wave functions 



F{V,6V) = \(*o(y)\* (y + 6V))\. 
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The fidelity metric g is a dimensionless, intensive quan- 
tity and can be defined as 
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where N is the number of sites. For finite-size systems, a 
decrease in the fidelity is a precursor to a quantum phase 
transition. This corresponds to a peak in the fidelity 
metric, which should diverge in the limit N — > oo for 
a second-order phase transition. For a transition to the 
CDW phase, this will coincide with a peak in 5(k = 0), 
which will similarly diverge. The fidelity metric also 
tracks level crossings well. At a level crossing, the wave 
functions |\&(V)) and \^(V + SV)) will be very different. 
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FIG. 3: (Color online) Total density of states and local density 
of states in the Haldane model for a site on the edge and in 
the bulk for (a) M = 0.0, (b) M = 0.5, (c) M = 1.0, (d) 
M = 1.1, (e) M = 1.2, and (f) M = 1.5. The calculations 
were performed on a N = 420-site cluster with ti — 1.0, 
£2 = 0.3, and <j> = t/4. The gap in the edge states opens 
at M c = 1.1. The vertical, gray line indicates the chemical 
potential. 



Because the overlap between the two states is negligible, 
the fidelity metric is sharply peaked. This, in conjunc- 
tion with a jump in the structure factor 5(k = 0), will 
distinguish a first-order phase transition. In the calcula- 
tions that follow, we take SV = 10~ 4 , which is sufficiently 
small to ensure results consistent with SV — > 0. 



III. NONINTERACTING SYSTEMS 

We first study the noninteracting fcrmions, which are 
exactly solvable, to compare our numerical results with 
the known solution. In addition, this limit allows for 
the study of large clusters and we can carefully address 
whether the small clusters that are accessible with the 
Lanczos algorithm are suitable for observation of a topo- 
logical insulator. 

As Haldane was able to show previously, 1 there is a 
phase transition between the topological insulator and 
the topologically trivial charge insulator as one increases 
the value of M. For M below the critical value of 
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FIG. 4: (Color online) Total density of states and local density 
of states in the Haldane model for a site on the edge and in 
the bulk for (a) M = 0.0, (b) M = 0.5, (c) M = 1.0, (d) M = 
1.1, (e) M - 1.2, and (f) M = 1.5. The calculations were 
performed on a N — 24-site cluster with the same parameters 
as in Fig. 3. The vertical, gray line indicates the chemical 
potential. Because of the finite size, the gap opens at smaller 
M and is, in general, larger than in Fig. 3. 
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FIG. 5: (Color online) Nearest-neighbor currents for the Hal- 
dane (left) and modified Haldane (right) Hamiltonians on a 
iV = 24-site cluster and parameters t\ — 1.0, t% = 0.3, and 
(j> = tt/4 with (a) M = 0, (b) M = 0.10, and (c) M = 1.5. In 
the modified Haldane model, a non-zero circulating current 
exists for M 7^ 0, i.e., the moment chiral symmetry is broken. 



M c = 3 3 / 2 t2sin0, the system is a topological insula- 
tor with an insulating bulk and a chiral metallic edge. 
As M approaches M c , the insulating gap in the bulk de- 
creases and closes at M = M c while the system becomes a 
semimetal. If one further increases M, an insulating gap 
develops in both the bulk and the edge. However, this 
insulator is topologically trivial and there are no gaplcss 
edge states with topological protection. 

Unlike an ordinary phase transition, this transition was 
beyond the framework of the Landau's theory of phase 
transitions. The two phases discussed above cannot be 
distinguished by any local parameter and no symmetry 
was broken across this transition. As a result, Landau's 
free energy cannot be defined. The only distinction be- 
tween these two phases are their topological property, 
which is described by a topological index known as the 
Chern number. 1,43 ' 44 Experimentally, one can distinguish 
these two phases by their edge properties. We also em- 
phasize here that this phase transition shows no discon- 
tinuity for any thermodynamical quantity and a similar 
example of a topological transition is the transition be- 



tween the strong and weak pairing in p x + ip y supercon- 
ductors (at least within the mean-field treatment). 44,45 

This topological transition can be observed by measur- 
ing the local density of states [Eq. (3)]. Since we deal here 
with finite-size systems, the energy spectrum is discrete. 
To compensate for this, we have broadened the delta 
function peaks with a Lorcntzian of width 7 = O.Olii. As 
shown in Fig. 3, for small M, the LDOS has no energy 
gap for a site on the edge, indicating clearly a metal- 
lic edge. However, for the bulk site, the same spectrum 
becomes gapped, reflecting an insulating bulk. 

To estimate the effects of finite-sized systems on this 
transition, we reduced the size of the cluster down to 
24 sites (12 unit cells). In Fig. 4, the bulk gap and the 
gaplcss edge states can still be observed, indicating that 
the topological property of the anomalous quantum Hall 
state is robust with regard to finite-sized clusters. A sim- 
ilar effect was also observed in a system fundamentally 
different from the model studied here. As reported in 
Ref. 46, the topological ordering in the Kitaev model 47 
was observed for small system sizes. While the topo- 
logical ordering found by Chen et al. has little direct 
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FIG. 6: (Color online) Total density of states and local density 
of states in the modified Haldane model for a site on the edge 
and in the bulk for (a) M = 0.0, (b) M = 0.5, (c) M = 1.0, 
(d) M = 1.1, (e) M = 1.2, and (f) M = 1.5. The calculations 
were performed on a TV = 420-site cluster with t\ = 1.0, 
ti — 0.3, and <f> = 7r/4. The gap in the edge states opens 
at M c = 1.1. The vertical, gray line indicates the chemical 
potential. 



relation with the results presented here, the fact that the 
same phenomenon is observed in these two systems sug- 
gests that robustness against finite-size effects might be 
a common property shared by different topological states 
of matter. 

We also computed the nearest-neighbor current 
[Eq. (6)] in the system, which is shown for M = 0, 0.1, 
and 1.5 in Fig. 5 (next-nearest-neighbor currents are not 
shown because they are an artifact of the complex hop- 
ping in the Hamiltonian and exist for all parameters). 
Note that the topological state has a clear chiral edge 
current (rotating in the clockwise direction), while the 
current in the bulk is negligible. However, in the topolog- 
ically trivial insulator with M > M c , the chiral current 
remains until M becomes large. This is because Hal- 
dane's model explicitly breaks chiral symmetry, which 
allows for a rotating edge current. Due to this effect, 
by measuring only the edge current, it is not sufficient 
to distinguish the topologically trivial insulator from the 
topological insulator. Indeed, the proper means of dis- 
tinguishing the topological state is via the LDOS. 



For comparison purposes, we also studied a modifi- 
cation to Haldane's model [Fig. 2(b)], where the two 
sublattices have opposite flux patterns. This model is 
topologically trivial in all of the parameter regimes. For 
M < M c , the system is a conductor (a Fermi liquid) at 
half filling but becomes a topologically trivial insulator 
for M > M c . At the transition point M = M c (which 
is identical to the transition point in Haldane's model), 
the system is a scmimetal. This transition was observed 
clearly in the LDOS (Fig. 6). In particular, the insulat- 
ing phase with M > M c was found to be gapped both 
in the bulk and on the edge, which, as expected, signifies 
a topologically trivial insulator. Moreover, the same be- 
havior can be found even for the small clusters with only 
24 sites (not shown). 

Although the modified Haldane's model is topologi- 
cally trivial, this system has its own significance in the 
study of the anomalous Hall effect. As in the Haldane's 
model described previously, this model also has no net 
magnetic flux. However, the flux pattern breaks the chi- 
ral and time- reversal symmetries for any M^0, even in 
the absence of the a net magnetic flux. These symme- 
try properties are in close analogy to the (non-quantized) 
Hall effect, although the magnetic field is zero after av- 
eraging over each unit cell in this case. Due to this sim- 
ilarity, the conducting phase in the modified Haldane's 
model with < M < M c is referred to as an anomalous 
Hall state. 43 ' 48 This anomalous Hall effect can be eas- 
ily observed by looking at the currents in the system (see 
panels on right side of Fig. 5). At M — 0, the chiral sym- 
metry is preserved, and the current on the left and right 
edges flows in opposite direction in order to preserve the 
chiral symmetry. However, for M > 0, the ground state 
of the system is the metallic anomalous Hall state, which 
breaks the chiral symmetry. Now, the currents at the 
opposite edge will not cancel each other. Instead, they 
flow in the same direction and form a rotating current 
circling around the system. We emphasize here that this 
edge current is not related with a nontrivial topological 
structure in the ground state wave function, and it only 
reflects the broken chiral symmetry. 43 ' 48 ' 49 



IV. INTERACTING SYSTEMS 
A. Spinless fermions 

Having shown that the topological insulator phase can 
be distinguished on small lattices, we now consider the 
properties of interacting systems. In this section, we con- 
sider spinless fermions with repulsive nearest-neighbor in- 
teractions (hard-core bosons are discussed in Sec. IV B). 
In the discussion that follows, we present detailed results 
for the parameters t\ = 1.0, t 2 = 0.3, <fr = 7r /4, and 
71/ = before discussing the results for general </>. 

At the non-interacting limit, we have shown above 
that, at M = 0, the system is a topological insulator 
and the two sublattices have the same occupation num- 
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FIG. 7: (Color online) Nearest-neighbor currents (left) for spinless fermions on a 24-site cluster with OBC, ii = 1.0, £2 = 0.3, 
4> — 7r/4, and interaction strengths (a) V = 0, (b) V — 1, (c) V = 2, (d) V = 3, and (e) V = 4. Local density of states (right) 
for the edge (OBC) and the bulk (shown for both OBC and PBC). The zeros of N e d se (u>) and A^uik, obc(^) ar e shifted upward 
for clarity, and the frequency to is shifted so that the chemical potential is at ui — (which is also indicated by the dashed 
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ber. In the strong-coupling limit V — > 00, however, the 
particles will all go to one of the two sublattices to avoid 
the energy penalty for occupying two neighboring sites. 
At half filling, one sublattice is fully filled while the other 
one is empty, resulting in a CDW insulating phase. Note 
with the formation of the CDW, translational symmetry 
remains but inversion symmetry was spontaneously bro- 
ken. Hence, we conclude that the system goes though a 
quantum phase transition as V is increased. More inter- 
estingly, the insulating phase in the weak-coupling limit 
is signified by a nontrivial Chern number. 1 ' 43,44 However, 
the CDW phase in the strong-coupling limit has a triv- 
ial Chern number (C = 0) and is topologically trivial. 
Therefore, the topology of the ground-state wave func- 
tion changes from nontrivial to trivial as V increases. 
This change in topology is referred as a topological tran- 
sition. Within the current knowledge of the topological 
insulators, it is not clear whether the quantum phase 



transition and the topological transition are, in general, 
two separate transitions or if they coincide to become a 
single phase transition. 

We first consider the topological transition by examin- 
ing the LDOS at the edge and in the bulk. In Fig. 7, we 
show local currents (left) and the local density of states 
(right) on a 24-site cluster with open boundary conditions 
for several interaction strengths. At V = [Fig. 7(a)], the 
system is a topological insulator, evidenced by the bulk 
insulating gap and the metallic local density of states 
on the edge. The circulating current on the edge is also 
consistent with the quantum Hall state. As the inter- 
action strength is increased, the gap in Abuik(w) closes 
while the states near the chemical potential in iV e dge(w) 
are depleted. For V = 2, 3, and 4 [Fig. 7(c)-7(e)], the 
edge states are fully gapped and the system is a triv- 
ial insulator. The local density of states in this phase 
is characterized by a near complete depletion of the bulk 




FIG. 8: (Color online) (a) Fidelity metric g as a function 
of interaction strength for various clusters with parameters 
ti = 1.0, ti = 0.3, tj> = tt/4, and SV = 10~ 4 . (b) Scaled 
CDW structure factor N *' 2v Scdw for the same parameters 
(shown for r\ = and yjv = 2). All clusters that contain the 
zone corner (k = K) as a valid k point exhibit a sharp, first- 
order transition at V c ~ 2.0. Without this k point, a smooth, 
continuous transition is observed. The inset shows the first 
derivative of the structure factor. The discontinuity (i.e., the 
peak) marks the transition point. 



states closest to the gap in the upper band, effectively re- 
sulting in a larger insulating gap in the bulk than on the 
edge. The topological transition takes place at Vr ~ 2. 
Finally, we note that the trivial insulator at intermediate 
coupling is characterized by small but non-zero edge cur- 
rents that travel in the opposite direction from the edge 
currents in the topological insulating phase. In the limit 
of V ^> ti, the edge currents vanish. 

Next, we consider the nature of the quantum phase 
transition which leads to the spontaneous breaking of 
inversion symmetry. Here, the momentum space prop- 
erties of the cluster are vital to the type of transition 
observed. The noninteracting band structure for this 
Hamiltonian features a Dirac cone at the zone corner 
K = (27r/3, 27r/v3). For an infinite system, this will 
always be a valid momentum point, whereas for small 
clusters this is not generally the case. Because the sur- 
face states in the band structure pass through this point 
in reciprocal space, the order of the transition has a clear 
dependence on the choice of clusters. This is reflected in 
the fidelity metric and the scaled CDW structure factor 
N i r / 2 "jS 'cdw i which are shown in Fig. 8 for clusters from 
N = 16 to N = 30 sites with periodic boundary condi- 
tions (see the Appendix for more details on the cluster 
shapes). If one studies clusters with a reciprocal space 
that does not contain the zone corner, then a broad max- 
imum in the fidelity metric and a crossing of the scaled 
structure factor is observed at V « 2.0. On just this 
information, one would conclude that there is a second- 
order phase transition at this interaction strength. How- 




FIG. 9: (Color online) Scaled structure factor N~~ ,/2v Scow 
for (a) 2 = 1 and y/v = 1.96 (rj = 0.04) and (b) Z = 2 and 
7/V = 2 (rj = 0). In both cases, the transition point is found 
to be V c = 2.1±0.1. In the insets, the horizontal axis is scaled 
as well, leading to the collapse of all data points into a single 
curve. 



ever, clusters that contain the K point (18A, 18C, 24C, 
24D, and 30A) exhibit markedly different behavior. In- 
stead, we find that the fidelity metric exhibits a very 
large and sharp peak (18C, 24C, and 30A) or discon- 
tinuity (18A and 24D) at V « 2.0. Furthermore, the 
CDW structure factor exhibits a jump at this interaction 
strength, indicating a first-order phase transition. In- 
deed, when we examine the first derivative of the struc- 
ture factor [see Fig. 8(b) inset], there is removable discon- 
tinuity at V c sa 2.0, which indicates that the topological 
transition and the quantum phase transition occur at the 
same interaction strength. 

For the clusters that do not contain the K point, 
the quantum phase transition is a conventional one (in 
contrast to the topological phase transitions discussed 
above). Because the inversion operator and the iden- 
tity operator form a Z2 group, the symmetry-breaking 
pattern of this quantum phase transition belongs to the 
Zi (also known as Ising) universality class. If the tran- 
sition is second-order, it is then expected to follow the 
scaling relation of the Ising model in d + z dimensions. 
Here d = 2 is the spatial dimension and z is the dynamic 
critical exponent, which is usually a positive integer indi- 
cating how many spatial dimensions the time dimension 
shall be counted as. 50 For z = 1, the critical point follows 
the three-dimensional (3D) Ising scaling, for z = 2 and 
above, the exponents are the mean-field exponents. 51,52 

We expect that the structure factor will scale according 
the following rule 



N^^Scuw = f[(V - K)iV 1/2 1 , 



(11) 



where N is the number of sites and 7 = v(2 — rf). For a 
3D Ising model, 51 (corresponding to the case of z = 1), 
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FIG. 10: (Color online) Nearest-neighbor currents (left) for hard-core bosons on a 24-site cluster with OBC, t\ = 1.0, ti = 0.3, 
cf> — 7r/4, and interaction strengths (a) V — 0, (b) V — 1, (c) 1/ = 2, (d) 1/ = 3, and (e) V = 4. Local density of states (right) 
for the edge (OBC) and the bulk (shown for both OBC and PBC). The zeros of Ned se (u>) and iVbuik, obc(^) are shifted upward 
for clarity, and the frequency uj is shifted so that the chemical potential (dashed vertical line) is at uj = 0. 



v = 0.7 and r\ — 0.04, while in four dimensions (4D) 
and above 52,53 (z > 2), v = 0.5 and r) = 0. As the 
proper value of z is still an open question, we show the 
rcscaled structure factor N~ j / 2v Scdw as a function of 
interaction strength for both of these cases in Fig. 9. 
If the transition is second order, the set of curves will 
converge at the quantum critical point V — V c . Here, 
the critical point obtained in both cases remains almost 
unchanged. Moreover, we show the universal scaling re- 
lation, Eq. (11), in the insets of Fig. 9 and observe the 
same qualitative behavior for both z = 1 and z = 2. As 
there is no discernible qualitative difference between the 
scaling with z = 1 and z = 2, it is clear that higher or- 
der corrections in the finite-size scaling are relevant. For 
the largest system size shown in Fig. 9, the linear system 
size is L — N 1 / 2 sa 4.9. Thus, significantly larger system 
sizes need to be studied in order to properly determine 
the value of z. 

In Fig. 1(a), we present a schematic phase diagram 



for this model in the V-(f> plane. Note that the phase 
boundary shown in the figure is determined for a 24-site 
cluster from the maximum in the derivative of S*cdw [see 
Fig. 8(b) inset]. At <f> = 0, there is no magnetic field in 
the system to break the time-reversal symmetry and, at 
weak coupling, the system is a semimetal, not a topo- 
logical insulator. Once time-reversal symmetry is broken 
(0 < < 7r), a topological insulator is observed for small 
V. There is a topological transition from the topological 
insulator, and this coincides with a first-order quantum 
phase transition to a topologically trivial Mott insula- 
tor. Furthermore, we note that the topological phase is 
strongest at <f> = it/2. 



B. Hard-core bosons 

Next, we examined the properties of the Haldanc 
model with hard-core bosons at half filling and nearest- 
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FIG. 11: (Color online) (a) Fidelity metric g as a function of interaction strength for hard-core bosons in various clusters with 
parameters t 2 = 0.3, <j> = tt/4, and SV = 1(T 4 . [(b)-(d)] Scaled CDW structure factor N" 1 ^" S C DVf for (b) 3D Ising (z = 1), 
(c) 3D XY (z = 1), and (d) 4D Ising / XY (z = 2). All of the curves intersect at the critical point V c = 3.27 ± 0.05 in (b) and 
(c) and V c = 3.28 ± 0.05 in (d). In the insets of (b)-(d), the horizontal axis is also rescaled and all of the data points collapse 
into a single curve. 



neighbor repulsive interactions, to check whether there 
are phases with non-trivial topological properties. We are 
motivated by the fact that in many instances hard-core 
bosons and spinless fcrmions are known to have similar 
properties. In one dimension, one can use the Jordan- 
Wigner transformation to map one onto the other, i.e., 
there is a one-to-one correspondence between them. In 
higher dimensions, one can also expect some similarities. 
For example, a simple band-structure calculation shows 
that, at half filling, the ground state of nonintcracting 
fcrmions on a hypercubic lattice in the presence of a 
staggered potential M, like the one in Eq. (1), is metal- 
lic for M = and becomes insulating for any nonzero 
value of M. For hard-core bosons, on the other hand, 
one has a superfluid (Bose-Einstein condensed) phase for 
M < M c and an insulator phase for M > M c , 54 ' 55 i.e., 
both models have equivalent ground-state phases. The 
only difference in this case is that for spinless fermions 
the metal-insulator transition always occurs at M c = 0, 
while for hard-core bosons the supcrfluid-insulator tran- 
sition occurs at a finite value of M c , which depends on 
the dimensionality of the system (M c — in one dimen- 
sion, as expected from the exact mapping). Hence, simi- 
lar physics as discussed in Sec. IV A for spinless fermions 
would not be unexpected in hard- core boson systems, 
where, of course, equivalent phases may just appear for 
different ranges of values of the Hamiltonian parameters. 

However, in Fig. 10, we show that this is not the case. 
For small interaction strengths, the system is a super- 
fluid with essentially no edge current. As V is increased, 



the system becomes an insulator simultaneously on the 
edge and in the bulk, indicating a direct transition to the 
topologically trivial Mott insulator. This Mott insulating 
state, like the spinless fermion model, exhibits a circulat- 
ing current on the edge and a depletion of the states 
in the bulk of the upper band at intermediate coupling 
strengths. Here the edge currents in the insulating state 
are larger than for spinless fermions and consequently 
vanish at larger interaction strengths. 

To determine the nature of the phase transition to the 
CDW phase we examined the fidelity metric and CDW 
structure factor for various system sizes. The peak in the 
fidelity metric [Fig. 11(a)] grows with the system size, 
which is an indicator of a second-order phase transition. 
Unlike the spinless fermion case in Sec. IV A, the uni- 
versality class of the transition is unclear. We have per- 
formed scaling analysis of the structure factor for both 
Ising and XY universality classes with dimension d + z. 
The critical exponents for Ising and XY scaling are sim- 
ilar in 3D and identical in 4D and above (the upper 
critical dimension of both models is 4). Because these 
exponents are similar and finite-size effects are present 
in our calculations, we cannot determine which univer- 
sality class, if any, is the appropriate one. To illus- 
trate this, we show the structure factor rescaled by a 
factor of AT-t/ 2 " in Figs. 11(b)- 11(d) for 3D Ising, 3D 
XY, 56 and 4D Ising/XY, respectively. There is no dis- 
cernible difference between the three cases. The curves in 
A r " 7 / 2ly 5cDW versus V for all three cases cross at roughly 
the same interaction strength (V c = 3.27±0.05 in 3D and 
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V c = 3.28 ± 0.05 in 4D) and the crossing point coincides 
with the peak in the fidelity metric. Moreover, when 
plotting N^/' 2u S C dw versus (V-V c )N 1 / 2u (see Fig. 11 
insets), all the resulting curves almost lie on top of each 
other. A more detailed study of this model will presented 
in future work. 36 



V. SUMMARY 

In this paper, we have presented the exact nonper- 
turbativc study of strong correlation effects in topolog- 
ical insulators. We showed that the clusters that can 
be studied with the Lanczos algorithm are sufficiently 
large to identify a topological insulator, clearly distin- 
guishing that the edge states are conducting and the bulk 
states are insulating. For spinless fcrmions in the Hal- 
dane model with repulsive nearest-neighbor interactions, 
we found that for weak interaction strengths the system 
is a topological insulator with circulating edge currents. 
For large interaction strengths, the system is a topologi- 
cally trivial CDW insulator. By using the fidelity metric 
and the CDW structure factor, we showed that the choice 
of the cluster is significant and that clusters must have 
a reciprocal space that incorporates the zone corner of 
the Brillouin zone as a valid k point. For these clusters, 
we have observed that the transition from the topologi- 
cal insulator to the topologically trivial CDW insulator 
is first order. 

We also investigated the properties of the Haldane 
model with hard-core bosons at half filling and <\> = 7r/4. 
Here we found that the bosons are in a supcrfluid phase 
for weak interaction strengths and a Mott insulator for 
large interaction strengths. The growth of the fidelity 
metric with system size is consistent with a second-order 
phase transition but scaling of the CDW structure factor 
cannot distinguish whether the universality class, if any, 
is Ising or XY. Additionally, no signature of topological 
order is present in this system. However, this model ex- 
hibits interesting critical behavior with a phase diagram 
that is currently being investigated. 36 

We note there that the interacting Haldane model by 
no means exhausts all possible types of interacting lat- 
tice Hamiltonians, which may show interesting interplay 
of strong correlations and topological textures, but it 
rather represents the simplest case where the asymptotic 
regimes are well understood. Therefore, it is conceivable 
that other models may host even more interesting tran- 
sitions and phases such as topological Mott insulators, 
which would represent the lattice analogs of fractional 
quantum Hall states. Some candidate states of this sort 
were previously discussed in Refs. 57-59, however no con- 
vincing evidence for their existence has been found in any 
model. We conclude by pointing out that our work has 
important consequences for the studies of topological in- 
sulating states in interacting lattice systems because it 
demonstrates explicitly that most hallmark features of 
topological insulators survive in small systems, which can 




FIG. 12: (Color online) Illustration of different clusters. Only 
the 18A, 18C, 24C, 24D, and 30A clusters contain the zone 
corner as a valid k point. The 24A cluster was studied with 
both open and periodic boundary conditions. The basis vec- 
tors for the lattice are shown in the bottom right of the figure. 



be analyzed using unbiased exact diagonalization meth- 
ods and could be realized in experiments with trapped 
ions or ultracold gases in optical lattices. 60 
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Appendix: Cluster Selection 

It is common in the Lanczos literature to diagonal- 
ize clusters with L x L sites in addition to other square 
clusters which completely cover the lattice. 61 For a hon- 
eycomb lattice, we generally follow a similar prescription 
of choosing a parallelogram or "tilted rectangle" to de- 
scribe the cluster. These vectors can be written in terms 
of the basis vectors ai and a 2 , shown in Fig. 12, that 
describe the underlying lattice 



Ai = miai +TO2a 2 , 
A 2 = niai + n 2 a 2 , 



(A.l) 



where mi, m 2 , m, and n 2 are integers. In Fig. 12, we 
explicitly show the shape of the clusters used in this work. 
Note that some of these clusters do not have all of the 
symmetry properties of the bulk. The lone exception to 
choosing a parallelogram to tile the lattice is the cluster 
designated 24A This hexagon-shaped cluster can tile the 
lattice and retains its complete symmetry. It is the only 
cluster studied with open boundary conditions. 
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As mentioned in the main text above, it is crucial that 
one distinguish between two different types of clusters, 
which depend on whether K — (2n/3, 2tt/^/S) is a valid 
point in reciprocal space. For this purpose, we provide a 
necessary and sufficient condition for the K point to be 
a valid momentum point, 



ax = 2mi/3 — m2/3 € Z, 
«2 = 2m/3 — ri2/3 S Z. 



(A.2) 



If both a.\ and 012 for a particular cluster are integers, 
then K is a valid point in the reciprocal space for that 
cluster. The number of sites in a cluster is N = 2|miri2 — 
m-2ji\ I and it can be easily checked this condition can 
only be satisfied when the number of sites are 6/, where 
I is a positive integer. For the clusters we studied, only 
clusters of N = 18, 24, and 30 sites can contain K point 
as a momentum point. Notice that a different choice of ai 
and a2 will change the condition described in Eq. (A.2). 
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